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ABSTRACT 

As the only directly imaged multiple planet system, HR 8799 provides a 
unique opportunity to study the physical properties of several planets in parallel. 
In this paper, we image all four of the HR 8799 planets at H-band and 3.3/im 
with the new LET adaptive optics system, PISCES, and LBTI/LMIRCam. Our 
images offer an unprecedented view of the system, allowing us to obtain H and 
3.3/im photometry of the innermost planet (for the first time) and put strong 
upper-limits on the presence of a hypothetical fifth companion. We find that all 
four planets are unexpectedly bright at 3.3yum compared to the equilibrium chem- 
istry models used for field brown dwarfs, which predict that planets should be 
faint at 3.3/im due to CH4 opacity. We attempt to model the planets with thick- 
cloudy, non-equilibrium chemistry atmospheres, but find that removing CH4 to 
fit the 3.3/im photometry increases the predicted L' (3.8/im) flux enough that 
it is inconsistent with observations. In an effort to fit the SED of the HR 8799 
planets, we construct mixtures of cloudy atmospheres, which are intended to rep- 
resent planets covered by clouds of varying opacity. In this scenario, regions with 
low opacity look hot and bright, while regions with high opacity look faint, sim- 
ilar to the patchy cloud structures on Jupiter and L/T transition brown-dwarfs. 
Our mixed cloud models reproduce all of the available data, but self-consistent 
models are still necessary to demonstrate their viability. 



1. Introduction 

Efforts are underway to characterize the first generation of directly imaged extra-solar 



plane ts. A principal focus has been the HR 8799 planetary system f JMarois et al.l 12008 



20101 ). which with four planets, is currently unique as a directly- imaged multiple planet 



^Thc LBT is an international collaboration among institutions in the United States, Italy and Germany. 
LBT Corporation partners are: The University of Arizona on behalf of the Arizona university system; 
Istituto Nazionale di Astrosica, Italy; LBT Beteiligungsgesellschaft, Germany, representing the Max-Planck 
Society, the Astrophysical Institute Potsdam, and Heidelberg University; The Ohio State University, and The 
Research Corporation, on behalf of The University of Notre Dame, University of Minnesota and University 
of Virginia. 
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system. Studying these planets simultaneously is particularly powerful given their connected 
formation histories and appearances. 

HR 8799 is a young, A5V star with a A Boo deficiency of heayy metals and 3 d i stinct 



circumstellar d ust structures (JMarois et al.l l2008t ICowley et al.l Il969l : iGray fc Kayd Il999 



Su et al.l I2OO9I ) . There is some disagreement about the age of the system. Traditional 
age-dating methods, such as galactic space motion a nd Hertzprung- Russell diagrana posi- 



tion suggest tha t HR 8799 has an age o f 20-160 Myr (IMoor et al.ll2006l : Marois et al.ll2008 



Hinz et al.ll2010l : IZuckerman et al.ll201ll ). while astroseismology estimates are more consis- 
te nt with ^1 Gyr , which would make the planets significantly more massive brown dwarfs 



( iMoya et al.ll2010l ). Interestingly, the d ynamical stability of the planets themselves places 
uppe r -limits on the masses of the pl anets (JFabrycky fc Murray-Clayll2010l : lMoro-Martin et al. 
2OIOI : ISudol &: Haghighipourll2012[ ). which directly con verts to a young system age, based 
on th e planets' photometry and evolutionary models (JBurrows et al.l 119971 : IChabrier et al. 
2000h . 



An important implication of the relative youth and low-masses of the HR 8799 planets 
is that their appearances and atmospheric properties might be different than field brown 
dwarfs, which can have the same effective temperatures as giant planets while being older 
and more massive. Brown dwarf spectra have been used as proxies for giant planet spectra 
to plan direct imaging surveys and to interpret early discoveries. However initial results 
show that there are several key differences between the atmospheres of giant exoplanets and 
brown dwarfs. 

For field brown-dwarfs, the L^T spectral type transit ion, occurs at ^1200-14 00 K, 



where dust clouds settl e/condense below the photosphere ([Saumon fc Marleyl l2008l ). and 
CO is converted to CH4 feurrows et"aDl2003[ ICeballe et al.ll2002[ l. For the HR8799 planets, 
clouds are suspended in t he photosphere at lower effective temperatures (900-1200 K) than i s 



typical for brown dwarfs ( ICurrie et al.kOllt iMadhusudhan et al.ll201ll : iBarman et al. 



2011al ). 



Additionally, there appears to be more CO than CH4 relative to equilibrium chemistry 
models, implying that convection is mixin g hot material into the photosphere fa ster than 
the C0f->CH4 reaction can re-equilibriate (JHinz et al.ll2010l : iBarman et al.ll2011al ). Similar, 



but more extreme results have been found for 2MASS 1207 b, a 5-7 M. 



Mjup T ^ Hya brown- dwarf primary ( IChauvin et al.ll2004j : ISkemer et al. 



J4 



companion to a 25 



2011 



Barman et al. 



201 Ibl ). Evidently the HR 8799 planets and 2MASS 1207 b look similar to L-type brown- 
dwarfs, despite having effective temperatures more consistent with T-type brown-dwarfs. 

Multi-wavelength photometry and spectroscopy are the keys to understanding the dif- 
ferences between brown dwarfs and giant planets. In particular, working over a broad wave- 
length range is critical for understanding clouds, chemistry, and the radiative budget of 
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extrasolar planets. The challenge of working over a broad wavelength range is that adaptive 
optics (AO) systems perform better at longer wavelengths where atmospheric turbulence is 
less severe. But background radiation increases at longer wavelengths, and most AO sys- 
tems have numero us warm optics, which can make working at long wavelengths impractical 



( lLloyd-Hartll2000l ) . The sweet-spot for most AO systems has typically been the near- infrared 



1-2.5/im), although for extrasola r planets, there is a benefit to working at L' (3.8/im) where 



the planet-star contrast improves (JHeinze et al.ll2008l ). 



The Large Binocular Telescope (LBT) AO system can increase the wavelength range 
over which we study extrasolar planets. With its 672-actuator deformable secondary mirror 
(installed on one side of the telescope at the time of our observations) the system produces un- 
precedented image quality and contrast at short wavelengths. And because it is a deformable 
secondary AO system, it has a minimal number of warm optics, so that background noise 
remains low at long wavelengths. 

In this paper, we present Large Binocular Telescope (LBT) First Light Adaptive Optics 
(FLAO) images of HR 8799 at H-band (L65/xm) and 3.3/im, detecting the four planets (b-e) 
at both wavelengths. Both images are superior to previous attempts at these wavelengths 
due to the high performance of the LBT's AO system. Our images are the first detections of 
HR 8799 e at H-band and 3.3/im and the first unambiguous detections of HR 8799 b and d 
at 3.3/im. In Section 2, we give a basic description of the instrumental setup for the FLAO 
system, the near-infrared imager, PISCES, and the mid-infrared imager, LMIRCam, which 
is a component of LBTI. Additionally, we describe our data reduction procedure, which is 
an implementation of the Locally Optimized Combination of Images algorithm (LOCI). In 
Section 3 we estimate photometry for the four planets, based on our LOCI reductions. In 
Section 4 we use our H-band image to search for additional companions interior to HR 8799 e, 
taking advantage of the unprecedented contrast afforded by the LBT AO system. In Section 
5 we present thick-cloud, non-equilibrium chemistry model atmospheres and mixed cloud 
atmospheres in an effort to explain the appearances of the HR 8799 planets. We conclude in 
Section 6 and make suggestions for futur e work characte r izing HR 8799 and other directly- 



imaged exoplanets. A companion paper, lEsposito et al.l ( l2012l ). describes the instrumental 



setup for the AO system and PISCES in detail, provides an independent analysis of the 
H-band data (along with new Ks-band data), and presents astrometry and a new orbital 
analysis of the system. 
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Observations 



2.1. PISCES H-band 



We observed HR 8799 at H-band (A =1.66/im; FWHM=0.29/im) with PISCES flMcCarthv et al 



20011 ) on UT Oct 16, 2011, during Science Verification Time for tfie LET First-Liglit Adaptive 
Optics System (FLAO). Tlie LBT's FLAO system (Pl-Simone Esposito) is a 672-actuator 
deformable secondary adaptive optics system tliat makes use of an innovative pyramid wave- 
front s ensor, producirig high Streh l-ratio, low background images over a broad wavelength 
range (JEsposito et al.l l2010l . 1201 if ). At the time of our observations, one adaptive optics 
system was installed on the right telescope, so for the observations presented in this paper, 
only one 8.4 meter mirror is used. PISCES (PI-Don McCarthy) is a 1-2.5/im imager with a 
Hawaii 1024x1024 HgCdTe array, which at the LET (single 8.25 meter aperture), critically 
samples a diffraction-limited point-spread-function (PSF) at H-band (with a plate scale of 
0.0193" per pixel). For the observations in this paper, PISCES was installed on the front 
bent- Gregorian focus of the LBT's right side to take advantage of the LBT's new AO system 
while facility infrared cameras are being delivered and commissioned. 

During our H-band observations of HR 8799, conditions were photometric and the nat- 
ural seeing, as measured by a DIMM (Differential Image Motion Monitor) on the telescope 
structure, was as good as ~0.9". We obtained 901 images with 2-second integrations, over 
the course of 2 hours (90 degr ees of sky- rotation) with the telescope rotator turned off (ADI- 
Angular Differential Imaging: iMarois et al.ll2006l ). PISCES' readout-time is 6 seconds, which 
means our observations were inefficient. However, high-contrast observations are usually lim- 
ited by PSF stability rather than photon-noise, so the long readouts had a negligible effect 
on our final results. The 2-second integrations saturated out to a radius of ~0.15". We were 
not able to obtain unsaturated H-band images of the star for photometry and astrometry 
with PISCES' shortest integration time, 0.8 seconds. 



Images were processed to remove cross-talk and persistence as described in lMcCarthy et al 



(I2OOII ) using the corquad correction software!. We then dark-subtracted, flat-fielded, and 
distortion-corrected the images. Finally, images were aligned by maximizing their cross- 
correlation. We processe d the aligned images u sing the LOCI algorithm (Locally Optimized 
Combination of Images; iLafreniere et al.l 120071 ) . which has been shown to produce higher 



contra st images than other algorithms, such as ADI. In the terminology of I Lafreniere et al. 
( 20071 ). we used A''^ = 1, Na = 300, fif = 1, and a 1-pixel subtraction-region jj The values for 



^http://aries. as. arizona.edu/^observer/dot. corquad. pisces 

■^The general idea of LOCI with ADI is as follows: The stellar PSF is removed in a set of subtraction 
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Na and g are adopted from lLafreniere et al.l (120071 ) and we use A^^ = 1 instead of A^5 = 0.5 to 



suppress self-subtraction, although we note that changing A^^ to 0.5 has a negligible effect on 
our photometry. Our implementation of the LOCI algorithm also includes a FWHM-sized 
mask around every subtraction-region, to further suppress self-subtraction. 

After reducing all 901 images with LOCI, we evaluated the noise of the images (standard 
deviation of pixel count) within an annulus from 0.3"-0.5" . The noise dramatically increased 
after the 500th image, corresponding to an increase in natural seeing. We reran LOCI using 
just the first 500 images (1000 seconds, 61 degrees sky-rotation), marginally improving our 
final image, which is shown in Figure [1] HR 8799 b, c, d and e are all clearly visible. 

The LOCI algorithm assumes accurate knowledge of the star's position. However, the 
stellar core saturated even with PISCES' fastest readout. We were able to estimate the 
stellar centroid to within 0.5 pixels (0.01") based on the circular symmetry of the PSF. We 
then reran LOCI with a set of different stellar centroid positions comprising a 1x1 pixel box 
and found no significant astrometric or photometric discrepancies between the results. 



2.2. LBTI/LMIRCam 3.3/xm 

We observed HR 8799 at 3.3/im (A =3.31/im; FWHM=0.40/im) with LBTI/LMIRCam 



(Large Binocular Tele scope Interferometer/L and M-band Infrared Camera; iHinz et al.ll2008 



Skrutskie et all l2010f ) on UT Nov 16, 2011. Although LBTI will eventually be used to 
combine the light from both LBT apertures, it was used in single-aperture mode for these 
observations, since only one adaptive optics system was operational. LBTI (Pl-Phil Hinz) 
consists of a beam combiner (UBC- Universal Beam Combiner), which combines the light 
from the two telescope tertiary mirrors, and a science camera (NIC-Nulling Infrared Camera), 
which itself contains a 2-5/im channel (LMIRCam) and an 8-13yum channel (NOMIC-NuUing 
Optimized Mid-Infrared Camera). LMIRCam uses a Hawaii 2RG 2048x2048 HgCdTe array, 
which oversamples A =2-5yum PSFs (with a plate scale of 0.0106" per pixel) when using only 
one telescope primary. 

Conditions during our 3.3/xm observations were photometric, with a natural seeing (as 



regions, which together, comprise the image. For each subtraction region, the stellar PSF is estimated 
by constructing a linear combination of individual exposures that minimizes noise within a corresponding 
optimization region, which is centered on the subtraction region, but is much larger. The optimization region 
is an annulus with an area of Na PSF-cores and a ratio between its radial extent and azimuthal extent of 
g. Images whose parallactic angle differ by less than an amount so that a source would move by Ng FWHM 
are excluded from the optimization. 
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measured by the telescope DIMM) of ~1.1". We obtained 1920 images with 3.8 second 
integrations over the course of 3.3 hours (110 degrees of sky rotation). The 3.8 second 
exposure saturated inside of 0.06". We also obtained short (0.15 second) unsaturated images 
for astrometry/photometry. About 1% of the images had bad tip/tilt residuals and were 
removed. 

Images were flat-fielded and globally-bad pixels (defined as pixels further than 15% from 
the median fiat) were replaced with the average of the 8-nearest good pixels. The images were 
then nod-subtracted with nods taken every 20 images. After these basic reduction steps, we 
found that some pixels had changing biases over shorter timescales than the nod-subtraction. 
We removed these with three steps: subtracting the 3(j-clipped median of each column (to 
remove column bias- level effects), median-binning the data with a 2x2 box (since the PSF 
is over-sampled by a factor of 3 for a single aperture telescope), and replacing isolated bad- 
pixels (4-sigma from their neighboring pixels) with the median of the neighboring pixels 
For the remainder of this paper, LMIRCam "pixels" refer to binned pixels, which have a 
plate-scale of 0.0212". 

Images were processed with LOCI, using the same parameters and implementation de- 
scribed in Section 12.11 We evaluated the noise inside a 0.2-0.4" annulus in each LOCI- 
processed image, and removed ~10% of the images, which had high noise, due to residual 
bad-pixels and/or sub-optimum AO performance. Our final image, which is shown in Figure 
m clearly shows planets b, c, d and e. 



3. Photometry 

Since LOCI self-subtracts, we calibrated our photometry by subtracting fake planets 
from the raw data at the positions of the detected planets, and rerunning the full LOCI 
pipeline. Best-fit photometry and error bars were determined by adjusting the fluxes of the 
fake planets to determine the range of values resulting in reasonable subtractions. For the 
PISCES H-band data, the star was saturated, so we calculated photometry for planets c, d 
and e with respect to HR 879 9 b, and conver t ed to absolute magnitudes by adopting the 



magnitude of HR 8799 b from iMetchev et al.l (120091 ). For the LMIRCam 3.3/im data, we 



calculated photometry for all four planets with respect to unsaturated images of the star 
obtained immediately after the saturated images used to detect the planets. We converted 



these to absolute photometry using HR 8799's absolute magnitude at 3.3/im from lHinz et al. 



(I2OIOI ). Errors on the LMIRCam absolute calibration are primarily the result of changing 
Strehl-ratios and telluric absorption variation throughout our observations. To test the mag- 
nitude of the Strehl-ratio variation, we did r = X/D aperture photometry on the unsaturated 



data and found a standard deviation of only ~2%. To test the magnitude of the telluric 
absorption variation, we compared the Airy pattern of the saturated and unsaturated images 
and found them to be consistent within ~5%. Combining these error terms, we adopt an 
absolute calibration error of 0.06 mags for the LMIRCam 3.3/im data. H-band and 3.3/im 
photometry of HR 8799 are presented in Table [TJ 



Our H- band photomet r y for H R 8799 c and d are consistent with the results of lMarois et al 



(120081 ) and Metchev et all (120091 ). and we detect HR 8799 e for the first time at H-band, 
finding it to be ~0.3 mags brighter than the next brightest plane t (HR 8799 c). Our 
3.3/xm data are som ewhat inconsistent with previous photometry from lHinz et al.l (120 lOf ) and 
Currie et al.l (l201l[ ). which are independent reductions of observa tions taken with M MT/Clio 
(note that the Clio and LMIRCam 3.3//m filter s are identica l) . iHinz et al.l (|2010[ ) reported 
detections of HR 8799 c and d but not b, and ICurrie et al.l (l201l[ ) reported detections of 



HR 8799 b and c but not d. The most substantial disparity is for HR 8799 b, for which 



Hinz et al.l ( l2010l ) reported an absolute magnitude upper-limit of 14.82, while ICurrie et al. 



(1201 ll ) reported a d e tectio n of 13.96±0.28. Our detection of 13.22±0.11 is closer to the re 
suit of ICurrie et al.l (12011) but is stil l brigh ter by a significant amount. We rea nalyzed the 



f inal r educed images from lHinz et al.l ( 120101 ) and find the upper-limit reported by lHinz et al. 



( 120101 ) to be erroneous (lik ely a typog r aphic error ). Our photo r netry for HR 8799 c is also 
brighter than observed by iHinz et al.l (l2010l) and ICurrie et al.l (l201l[ ) and our photometry 
of HR 8799 d is brighter than observed by iHinz et al.l ( 120101 ). We present a comparison of 
the MMT/Clio image and our new LBT/LMIRCam im age in Figure El The fact that our 



photometry is somewhat inconsistent with the values of iHinz et al.l (|2010[ ) and ICurrie et al. 



( 2011 ) can likely be explained by the non-photometric conditions reported in iHinz et al. 



(120101 ). varying AO performance (caused by the non-photometric conditions) and overly- 
optimistic systematic and/or measurement error analyses considering the quality of the Clio 
data. Variability is unlikely to be a factor, given that it has not been reported in any other 
photometric band and the amplitude of variability needed to rectify the disparity is quite 
large. 



4. Constraints on a Hypothetical Fifth Planet, HR 8799 f 

While HR 8799 is known to have four giant planets at wide separations, the inner 
system might have one or more companions that have not yet been discovered. Additional 
companions would chall enge formation models, which already have a hard time explaining 
the outer four planets (JDodson- Robinson et al.l 120091 : iKratter et al.ll2010l ). A fifth planet 
would also complicate dynamical stability analyses, which require mean-motion resonances 
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and planet masses that are on t he low end of the range predicted by evolutionary models 



and independent age estima t es (IFabrycky fc Murray-Clayl l20ld : iMoro-Martm et al.l 12010 



Sudol fc Haghighipourll2012l ). iHinkley et al.l ( 20111 ) used non- redundant masking to rule out 
the presence of massive inner companions from ~0.01"-0.5", but only for objects significantly 
more massive than the four known planets. Here we use our H-band image to search for 
planetary mass companions interior to HR 8799 e. 

We evaluate our ability to detect a close-in companion by making a contrast curve of our 
residual H-band image (after the 4 planets have been removed) jj We smooth the image with 
a Gaussian that is the same size as our diffraction-limited PSF, and calculate the standard 
deviation in 1-pixel annuli. Counts are converted to photometry using the peak-fiux (central 
pixel) of planet 'b' from the smoothed image. We also correct for self-subtraction, which 
is measured by inserting fake planets into the raw data at various radii. Our 5a contrast 
curve is shown in Figure |3l HR 8799 b-e are shown as diamonds. The vertical dashed line 
denotes the separation of the 2:1 mean-motion resonance with HR 8799 e (assuming a face-on 
circular orbit for simplicity). If there is a massive interior planet, it is likely to be in a stable 
resonance, as has been found for pairs of the outer four planets. We find no fifth planet 
at or exterior to the 2:1 resonance with HR 8799 e, with limits down to the approximate 
brightnesses (and by extension, masses) of the inner three planets. As a check, we insert a 
fake 'e'-like planet into our raw-images at a separation of 0.235" , the approximate position of 
the 2:1 resonance with 'e'. Our reduced image (FigureH]) shows that we would have detected 
"planet f" anywhere exterior to the 2:1 resonance, and that there are no residuals brighter 
than 'f ' in the image. Note that we do not repeat this analysis with the 3.3/xm image because 
it is not as sensitive to additional companions as the H-band image. 



5. Multiwavelength Modeling of the HR 8799 planets 



Multiwavelength photometry and spectroscopy have improved our understanding of 
the physical properties of the HR 8799 planets, which are not well fit by th e same mod- 



els th at have been used to interpret the properties of field brown dwarfs. iMarois et al. 



(120081 ) first noted that the HR 8799 planet SEDs were best fit by model atmospheres with 
T=1400-1700 K, while their luminosities were more consistent with T=800-1000 K. The 



''Note that we use an H-band image construct ed with Ns = 0.5 instea d of A^5 = 1.0 (which was described 
in Section [3]) to maximize S/N, as described in iLafreniere et al.l (|2007l) . The Ns = 0.5 contrast curve is 



~0. 1-0.2 mags better than the Ns — 1.0 contrast curve from ~0.2-0.3", and then becomes progressively 
worse (up to ^^0.5 mags) outside of 0.3". For this section, the inner regions are more important, so we use 
the Ns = 0.5 contrast curve. 



-10- 



disparity was driven by the planets' faint/red appearance with respect to models, and a 
lack of methane absorption in narrow-band 1.59/1.68 CH4 photometr y. Subsequent 3.88- 
4.10yum spectr oscopy of HR 8799 c was inconsistent with both COND (JBaraffe et al.ll2003[ ) 
and DUSTY (jChabrier et al.l 120001 ) atm ospheric models, w hich the authors interpreted as 
evidence for non-equilibrium chemistry (jjanson et al.l 120101 ). 3-color photometry in the L 
and M bands were also inconsistent with equilibrium chemistry atm ospheric models, and 
in particular showed a relative lack of methane absorption at 3.3/im (IHinz et al.l 120101). 
and K spectroscopy o f HR 8799 b also s howe d a methane deficiency (iBowler et al.l 12010 



Barman et al.l l2011al ) . ICurrie et al.l ( 120111 ) and iMadhusudhan et al.l (1201 ll ) were able to pa- 
rameterize thick cloud models to fit multiwavelength photometry for HR 8799 bed. However, 
th e models were not able to reproduce the 3.3/im photometry or the subsequent spectroscopy 
of iBarman et al.l (l2011al ). who fit both the photometry and spectroscopy of HR 8799 b with 
models that incorporated clouds and non-equilibrium chemistry. 

The combined result of these studies is that HR 8799 b, c and d have non-equilibrium 
C0^CH4 chemistry and cloudy atmospheres at low effective temperatures where, in field 
brown dwarfs, the clouds are thought to have settled below the photosphere. Two color- 
magnitude diagrams, shown in Figure El demonstrate these effects. On the left, is an H- 
K versus H color-magnitude diagram, which shows the M— )-L— >T sequence of field brown 
dwarfs. L dwarfs are characterized by cloudy atmospheres, while T dwarfs are characterized 
by cloud-free atmospheres. The intermediate region is the L^T transition, where objects are 
thought to have patchy clouds or clouds that have partially descended below the photosphere. 
The HR 8799 planets appear to be an extension of the field L-dwarf sequence, i.e. they have 
clouds at an effective temperature where the field brown dwarfs are transitioning to cloudless. 
The right-hand side of Figure [5] shows a 3.3/im-L' versus L' color magni tude diagram, which 
includes a sequence of chemical equilibrium, thick cloud models from IMadhusudhan et al. 
(120111 ). The HR 8799 planets are all much brighter at 3.3/im then predicted by the models, 
implying a relative absence of CH4, which is a strong absorber at 3.3/im. 

In the following sections, we present model atmospheres of the HR 8799 planets in an 
attempt to reproduce their photometry. A summary of the available photometry is presented 
in Table 121 colors are presented in Table [31 and a listing of all the models used in this paper is 
presented in Table HI For all mod el comparisons in t his paper, we convolve the model planet 
atmosphere and a model of Vega (ICohen et al.lll992l ) with filter profiles to produce predicted 
magnitudes in the different filters, which are then compared to the measured photometry. 
The filter profiles have all been multiplied by a model telluric atmospherCl. For most filters. 



^http://www.gemini.edu/sciops/ObsProcess/obsConstraints/atm-models/cptrans^m_43_15.dat 
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this step has a neghgible effect, because the telluric atmosphere has a ffat transmission 
profile, but Earth's atmosphere has a large transmission slope in the 3.3/im filter which 
changes the effective wavelength of the filter, and in turn, changes the model magnitudes by 
~0.1 mags. For the sake of plotting the model fits, we convert the filter magnitudes to Jy 
by convolving the Vega model with filter profiles to produce zero-point fluxes. For HR 8799 
b, we also include the H and K spectroscopy from iBarman et al.l (J2011al ). The HK spectrum 
comparison is made by smoothing the planet model atmosphere by the published spectrum's 
resolution (0.01/im) and directly comparing to the observed spectrum. The absolute fluxes 
of the H and K spectroscopy are tied to the H and K photometry, so for our comparison, we 
allow the overall brightness of the two spectra to vary, and only fit the shapes of the spectra. 



5.1. HR 8799 b 



HR 8799 b is the most challenging atmosphere to explain because it is the coolest of the 
four planets and is thus the largest outlier in the color-magnitude diagrams shown in Figure 
m Additionally, the HR 8799 b atr nosphere is the most constrained of the four planets 
due to the H and K spectroscopy of iBarman et al.l ( l2011al ). Our new 3.3/im photometry 
of HR 8799 b is s ignificantly brighter (by ^100%) than the pr e yious published value by 



Currie et al.l (120 ll[ ). As described in Section [31 the ICurrie et al.l ( 120111 ) 3.3/im photometry 



is likely erroneous, due to non-photometric conditions and/or overly optimistic error bars. 
There is the possibility that HR 8799 b is extremely variable at 3.3/im, but we consider that 
scenario unhkely due to the many times HR 8799 has been observed at other wavelengths 
where there has been no evidence of such large variability. 



Before proceeding with our modeling , we examine IBarman et al.l ( l2011af )'s recent hy- 
pothesis regarding the SED of HR 8799 b. IBarman et al.l (2011a) were able to reproduce all 
existing photometry and spectroscopy of HR 8799 b by using a non-equilibrium C0^CH4 
cloudy atmosphere with Teff=1100K and log{g)=3.5. However, their model is inconsis- 



tent wi th interior evolutionary models, which predict a larger object radius. IBarman et al. 
(l2011al ) also present model atmospheres that obey the predictions of the evolutionary tracks, 
with a best-fit model that has T=896 K and Z=1.0 metallicity. This model demonstrates 
that it might be possible to explain HR 8799 b's appearance with a combination of clouds, 
non-equilibrium chemistry and higher than solar metallicity. However, the model does not 
fit all of the existing photometry, and our new 3.3/im photometry makes the fit significantly 
worse. 
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5.1.1. Non-equilibrium CO^CH^ Chemistry 



We start our modeling by using the parameterized thick-cloud atmospheres of lMadhusudhan et al 



( 120111 ) and adjusting the CO and CH4 mixing ratios. These models assume evolutionary 
track radii, but make no atte mpt to self-consistently explain the CO and CH4 mixing ratios, 
which iBarman et al.l ( l2011al ) produce with turbu lent mixing. In Figu r e El w e plot the best 
fit chemical equilibrium model of HR 8799 b from iMadhusudhan et al.l (120111 ). as well as two 
models that have suppressed CH4 and enhanced CO with respect to the equilibrium models 
(by factors of 10 and 100 respectively). Suppressing CH4 and enhancing CO has a negligible 
effect on the near-infrared (zJHK) photometry, but it greatly affects the HK spectroscopy, 
favoring the 100 xCO, 0.01 XCH4 model. Non-equilibrium chemistry dramatically affects the 
3-5/im SED, where a lack of CH4 makes the object brighter in the 3.3/im and L' filters and 
excess CO makes the object fainter in the M-band filter. The 100 x CO, 0.01 XCH4 model 
fits the 3.3/im photometry, and comes closest to fitting the M-band photometry (further 
increasing the CO mixing ratio would improve this fit). However, the lack of CH4 strongly 
increases the predicted flux in the L' filter, making it incompatible with HR 8799 b's observed 
L' photometry. Based on this analysis, it seems unlikely that a cool (850 K) atmosphere 
with non-equilibrium chemistry could explain HR 8799 b's 3.3/im-L' color. However, there 
remains the possibility that more complex radial profiles of non-equilibrium chemistry and 
clouds, could explain all of the data. 



5.1.2. Mixed Cloud Atmospheres 



Based on Figure |5l HR 8799 b has colors reflective of a ~1300 K atmosphere, but 
a luminosity consistent with an ~850 K atmosphere. In lieu of resolving the difference 
by assuming a small (unphysical) radius, we consider the possibility that the planet emits 
non-isotropically, with bright and dark sections, such that the bright sections dominate the 
shape of the SEP . Bright and dark regions have been observed on Jupiter in the mid- infrared 
( IWestphallll969[ ) . and several s tudies seeking to expl ain the L^T transition have used hybrid 
cloudy/cloud- free models (e.g. iBurgasser et al.ll2002l ). For our purposes, mixtures of standard 
cloudy and cloud-free models are unlikely to explain the appearance of HR 8799 b because 
its H-K color is redder than the L-dwarf sequence implying that it is more cloudy, not less 
cloudy, than the other L-dwarfs. 

An alternative is that HR 8799 b has a mixture of clouds, such that the whole planet 
is cloudy, but with regions that have thicker clouds where the planet appears darker. A 
coarse approximation of this phenomenon i s to linearly combine c loudy models of different 
effective temperatures (analogously to how IBurgasser et al.l (120021 ) and others have linearly 
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combined cloudy and cloud-free models of different effective temperatures). This method is 
somewhat non-p hysical due to the fa ct that the two atmospheres have different temperature- 
pressure profiles ( iMarley et al.ll2010[ ). Linearly combining models with a shared temperature- 
pressure profile but different cloud structures would be more correct, but is beyond the 
scope of this paper. In Figure [TJ we present an example of a hybrid atmosphere that is 
93% Teff=700 K, 'A'-ty pe cloudy and 7% r^ fF =140 K, 'AE'-type cloudy (A and AE cloud 
pr ofiles are descri b ed in iMadhusudhan et al.l (120111 ) and further model details are presented 



m 



Burrows et al.l ( l2006l )). The hybrid model adequately fits all of the photometry except 



for M-band, which can be explained with further increased CO absorption. The HK model 
spectrum is muted compared to the data, but its bulk shape is generally correct (i.e. it does 
not show strong CH4 absorption, as would be expected for a cool object). We note that 
this model is meant to be representative, but that a true mixed-cloud atmosphere should be 
calculated self-consistently. 



5.2. HR 8799 c, d and e 



We repeat our analysis of HR 8799 b (Section 15.11) for HR 8799 c and d, again starting 
with the best-fit, thick-cloud models from IMadhusudhan et al.l (120 ll[ ). A comparison of non- 
equilibrium C0-H-CH4 models is shown in Figure |H] and example mixed cloud atmospheres 
are shown in Figure IHi Our conclusions for HR 8799 c and d are similar to our conclusions 
for HR 8799 b: we are unable to fit the 3.3/xm-L' colors of HR 8799 c and d with cloudy/non- 
equilibrium chemistry models, and mixed cloud atmospheres do a reasonable job fitting all 
of the data. We purposely construct the mixed cloud atmospheres from the same two model 
atmospheres used to make the mixture for HR 8799 b (see Section [5. 1.2p . In this scenario, 
we find that HR 8799 c and d have higher fractions of the warm atmosphere than HR 8799 
b, explaining their higher luminosity. 

HR 8799 e has less data than the outer 3 planets, but the existing data is consistent 
with the photometry of HR 8799 c and d. With the addition of our new H-band and 3.3/im 
photometry, HR 8799 e has now been studied at four wavelengths, and we can proceed with 
atmospheric modeling, based on lessons learned from HR 8799 b, c and d. We begin by using 
the thick-cloud models from IMadhusudhan et al.l ( 1201 ll ) to look for a good fit of the H, K 
and L' photometry. We find that a 1000 K, log(g)=4.0 model fits well (shown in Figure [TU]) . 
Based on the small number of data points, degeneracies between different cloud properties, 
surface gravity and effective temperature are not explored, but it is reasonable to assume the 
cloud properties and surface gravity will be similar to the other HR 8799 planets. Figures 
ITT] and [12] show non-equilibrium chemistry models and mixed cloud models. As was true for 
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the outer three planets, we are unable to fit the available photometry with non-equilibrium 
chemistry models, and a mixed-cloud model can be made to fit. 



5.3. The HR 8799 Planets in Aggregate 

The HR 8799 system provides a unique laboratory for simultaneously studying mul- 
tiple coeval planets. In the previous sections, we have modeled the planets individually. 
Comparing the four planets provides additional insight. 

HR 8799 c, d and e are brighter than HR 8799 b in all published photometry, but 
have almost exactly the same colors as 'b' throughout their measured SEDs (see Table [3l 
the colors for all four planets are the same, within errors, except for J-H for HR 8799 c, 
where J is abnormally bright). Give n that the planets ar e coeval, they are expected to all 



have similar radii (to within ~10%; iBurrows et al.l 119971 ). Therefore, the fact that 'c', 'd 



and 'e' are brighter than 'b' implies that they have higher effective temperatures (based on 
L = An R'^ aT^r r) . However, the fact that all of the HR 8799 planets have the same colors 
suggests that the physical properties of their atmospheres are similar, despite their different 
effective temperatures. In field brown dwarfs, objects with different effective temperatures 
(over the range probed by the HR 8799 planets) have different physical properties and 
different SED colors. For the HR 8799 planets, this does not appear to be the case|j 

The fact that 'c', 'd' and 'e' have almost the exact same colors as 'b' is circumstantial 
evidence for the mixed-cloud atmospheres. In the mixed-cloud scenario, HR 8799 b has a 
lower effective temperature than the other HR 8799 planets, but a similar physical temper- 
ature in the bright, emitting regions of its atmosphere. In Figures [71 M and [T21 we have 
have shown that all of the HR 8799 planets can be fit by a mixture of 1400 K and 700 K 
cloudy atmospheres (which we have been using as an approximation for two different cloud 
structures, one of which has much lower opacity). The difference between the cooler planet, 
'b', and the hotter planets, 'c', 'd' and 'e', is the mixing ratio of the 1400 K and 700 K 
models, i.e. the covering fraction of the different cloud types. 



The non-equilibrium chemistry models (with reasonable radii) of iBarman et al.l ( l2011af ) 



fit the HK spectroscopy better than our makeshift mixed-cloud models, bu t they do not 



reprodu ce the broad-band colors as well, in particular from 3.3-3.8/im. The iBarman et al. 



f l2011af ) models do provide physical motivations for non-equilibrium chemistry (turbulent 



an exar 



mple, iDupuy fc Liul ( 2012 1) find AK-L'/AK=0.4 over a large effective temperature range, 



whereas HR 8799 b has a AK-L'/AK=0.03 with respect to the average values of HR 8799 c, d and e 
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radial mixing), while our mixed cloud models are based on analogy with physically quite 
different systems (Jupiter and L/T transition brown-dwarfs). Self-consistent modeling is 
still necessary to determine if mixed cloud atmospheres are a viable explanation for the HR 
8799 planets. In any case, it appears detailed modeling of increasingly complex cloud physics 
and chemistry will be necessary to explain the true nature of the HR 8799 planets. 



6. Summary and Conclusions 

We have directly imaged the HR 8799 planetary system, detecting all four planets at 
H-band and 3.3/im with the LBT's First Light Adaptive Optics system. The images are of 
unprecedented quality allowing us to rule out the presence of a massive (HR 8799 cde-like) 
planet exterior to HR 8799 e's 2:1 inner resonance (H-band 5-cr contrast of 11.6 magnitudes at 
0.235"). We detect HR 8799 e at H-band, for the first time, and find that it is approxirn ately 



as bright as HR 8799 c and d. Combined with Ks and L' data (JMarois et al.l l2010l ) and 
our new 3.3/im data, this indicates that HR 8799 e has similar atmospheric properties to 
HR 8799 c and d. The planets are all brighter than expected at 3.3yum, where equilibrium 
chemistry models predict CH4 opacity should make the planets faint. 



We model the HR 8799 planets with thick-cloudy atmospheres (JMadhusudhan et al. 



20111 ) and allow the CO and CH4 mixing ratios to vary arbitrarily. The models that fit 
our 3.3/im data (which have very little CH4) substantially over-predict the planets' fluxes 
at L' (3.8/im). Hotter atmosphere s (>1300 K) have a sim ilar 3.3/im-3.8/im color as the HR 



8799 planets, and the shape of the iBarman et al.l (J2011al ) HK spectrum is also well fit by an 



atmosphere that is significantly hotter than indicated by the HR 8799 planets' bolometric 
magnitudes. As a result, we consider the possibility that small sections of the planets' 
atmospheres are hot (>1300 K), dominating the shape of the SEDs, while the majority 
of the planets' atmospheres are cooler and do not produce much flux. The temperature- 
pressure profile must be the same between the "hot" and "cool" regions, so the physical 
difference would be that the "cool" regions have increased cloud opacity. Since the SED is 
consistent with a cloudy atmosphere, the "hot" regions must also be cloudy, so the combined 
atmosphere is comprised of mixed clouds, some of which are thicker than others. Our mixed 
cloud models are able to fit all of the HR 8799 data. However, we caution that our models 
are not fully self-consistent and that more theoretical work is necessary to validate our 
hypothesis. 

The HR 8799 planets have unusual SEDs that are not well-fit by the same models that 
have been used to fit field brown dwarfs. From an observational standpoint, HK spectroscopy 
and 3.3/im-L' colors have been particularly powerful in ruling out model atmospheres. HK 
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spectroscopy of HR 8799 cde, and other directly-imaged planets in general, is critical to 
our understanding of clouds and non-equilibrium chemistry. We also note that given the 
wide range of 3-4/im SEDs predicted by the models in this paper, it would be very useful to 
obtain low-resolution spectroscopy of the HR 8799 planets in this range. From a theoretical 
standpoint, our new 3.3/im photometry is a challenge even for non-equilibrium chemistry 
models, which predict bright 3.3/im photometry. Mixed cloud models are one possible way 
to flatten out 3.3/im-3.8/im photometry and hide CH4 opacity. 

The authors thank Piero Salinari for his insight, leadership and persistence which made 
the development of the LET adaptive secondaries possible. We also thank the many dedi- 
cated individuals who have worked on LBTI, LMIRCam, PISCES, and the AO system over 
the years. The Large Binocular Telescope Interferometer is funded by the National Aero- 
nautics and Space Administration as part of its Exoplanet Exploration program. LMIRCam 
is funded by the National Science Foundation through grant NSF AST-0705296. 
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Table 1. LBT Photometry of the HR 8799 Planets 



Planet AH Mag with HR 8799 b Absolute H Mag AS-S^m Mag with HR 8799 Absolute 3.3/im Mag 

HR 8799 b 15.08±0.13'' 10.97±0.10 13.22±0.11 

HR 8799 c -0.90±0.05 14.18±0.14 9.97±0.10 12.22±0.11 

HR 8799 d -0.85±0.2 14.23±0.2 9.77±0.10 12.02±0.11 

HR 8799 e -1.2±0.2 13.88±0.2 9.87±0.20 12.02±0.21 

''HR 8799 b absolute H-band photometry from lMetchev et alj 1 120091) . Absolute H-band photometry for the other planets 
is with respect to HR 8799 b. 



Table 2. Photometry of the HR 8799 Planets 



Planet 



z 
(1.03Atm) 


J 

(1.25/im) 


H 

(1.63/im) 


CH4S 
(1.59/im) 


CH4I 
(1.68/im) 


Ks 

(2.15/im) 


3.3 

{3.3^m) 


L' 
{3.8^lm) 


M 
{4.7/im) 


18.24±0.29 


16.30±0.16 
14.65±0.17 
15.26±0.43 


15.08±0.13 
14.18±0.14 
14.23±0.2 
13.88±0.2 


15.18±0.17 
14.25±0.19 
14.03±0.30 


14.89±0.18 
13.90±0.19 
14.57±0.23 


14.05±0.08 
13.13±0.08 
13.11±0.12 
12.93±0.22 


13.2±0.11 
12.2±0.11 
12.0±0.11 
12.1±0.21 


12.66±0.11 
11.74±0.09 
11.56±0.16 
11.61±0.12 


13.07±0.30 
12.05±0.14 
11.67±0.35 



HR 8799 b 
HR 8799 c 
HR 8799 d 
HR 8799 e 
reference (b,c,d,e) 



3,3,3 



4,1,1,1 



3,3,3 



3,3,3 



3,3,3,5 



1,1,1,1 



3,3,3,5 



6,6,6 



Ref erences. — (1) This work; (2) ICurrie et al.l l l201lh : (3) iMarois et al.N200Sh ; f4') lMetchev et al.l ll2009h ; (5) iMarois et al.N20ld) : (6) iGalicher et al.l 
l|201ll ') 



00 
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Table 3. Colors of the HR 8799 Planets 



Planet J-H H-Ks Ks-L' 3.3Atm-L' L'-M 

HR 8799 b 1.22±0.20 1.03±0.15 1.39±0.13 0.54±0.15 -0.41±0.31 

HR 8799 c 0.47±0.21 1.05±0.16 1.39±0.12 0.46±0.14 -0.31±0.16 

HR8799d 1.03±0.46 1.12±0.23 1.55±0.20 0.44±0.19 -0.11±0.37 

HR 8799 e 0.95±0.29 1.32±0.25 0.49±0.24 
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Table 4. Atmospheric Models Used in This Paper 



Figure # 


Te// (K) 


log(g) 


Cloud Type 


Chemistry 


m 


850 


4.3 


AE.60 


equilibrium 
10 xCO, O.IXCH4 
lOOxCO, O.OIXCH4 


El 


1400 
700 


4.0 


AE.60 
A.lOO 


equilibrium 



H 1000 4.2 AE.60 equilibrium 

10 xCO, O.IXCH4 
lOOxCO, O.OIXCH4 
900 3.8 AE.60 equilibrium 

10 xCO, O.IXCH4 
lOOxCO, O.OIXCH4 



m 


1400 
700 


4.0 


AE.60 
A.lOO 


equilibrium 


\M 


1000 


4.0 


AE.60 


equilibrium 


\m 


1000 


4.0 


AE.60 


equilibrium 
10 xCO, O.IXCH4 
lOOxCO, O.OIXCH4 


m 


1400 
700 


4.0 


AE.60 
A.lOO 


equilibrium 



Note. — Cloud types ('A' and 'AE) refer to cloud thickness and the associ- 
ated numbe rs ('60' and ' 100') refer to modal dust grain size, as parameterized 
in Madhus udhan et al.l 11201 j) . Multiples of CO and CH4 are with respect to 
equilibrium chemistry model s. All models assume solar metallicity. Further 
model details are discussed in lBurrows et al.l 1120061 ) and lMadhusudhan et al.l 
1I2OIII') . 
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Fig. 1. — LBT First Light AO images of the HR 8799 planetary system at H-band and 
3.3/im. These images comprise the first detection of HR 8799 e at either wavelength, and 
the first unambiguous detections of HR 8799 b and d at S.Sfim. 
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Fig. 2. — A comparison of the MMTAO/Clio 3.3/im image from lHinz et al.l (J2010l ) and our 
new LBTAO/LMI RCam 3.3/im imag e. Our 3.3/im photom etry is somewhat inconsistent 
with the findings of iHinz et al.l ( l2010l ) and lCurrie et al.l ( 1201 ll ). who separately analyzed the 
MMT data. Based on the relative quali ty of the iraages, it is l ikely that the dispa rity is the 
result of overly-optimistic error bars by lHinz et al.l (l2010[ l and lCurrie et al.l (l201l[ l. 
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Fig. 3.— H-band contrast curve for the LBTAO/PISCES H-band image of HR 8799, with 
the four planets shown. Also shown is the position of the 2:1 orbital resonance with HR 
8799 e (assuming, for simplicity, a face-on, non-eccentric orbit). If there is a massive inner 
planet, it is likely to be in a stable resonance, as has been found for the outer companions. 
Based on the contrast curve, we would have been able to detect a planet at HR 8799 e's 2:1 
inner resonance, if it were approximately as bright (massive) as HR 8799 cde. Note that the 
contrast curve shows a d ark hole inside of ^ 0.6", which is a predicted feature of high-order 
adaptive optics systems ( iMalbet et al.lll995[ ). 
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Fig. 4. — Image of the HR 8799 system with a fake planet, 'f ', added in at the approximate 
location of HR 8799 e's 2:1 orbital resonance. The fake planet, which is the same brightness 
as HR 8799 e, was added into our individual raw frames and is easily recovered by our LOCI 
pipeline. There are no residual point-sources as bright as the fake planet at or exterior to 
its position at the HR 8799 e 2:1 orbital resonance. 
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Fig. 5. — LEFT: H vs. H-K color -magnitude diagram showing the M— >-L ^T spectral- type 
transition for field brown dwarfs ( jLeggett et al.ll2002l : iKnapp et al.ll2004j ) and the HR 8799 
planets. The HR 8799 planets appear to be an extension of the L-dwarf sequence, implying 
that they have cloudy atmospheres at lower effective temperatures than are typical for cloudy 
field brown dwarfs. 

RIGHT: L' vs. 3.3/xm- L' color-magnitude diagram showing equilibrium chemistry, thick- 
cloud atmospheres from lMadhusudhan et al.l (120111 ) and the HR 8799 planets . The HR 8799 
planets are all brighter at 3.3/xm than predicted by the iMadhusudhan et al.l (120111 ) models, 
implying a lack of CH4, which is a strong absorber at 3.3/im in the equilibrium chemistry 
model atmospheres. 
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Fig. 6. — Photomet ry and spectroscopy of HR 8799 b (red error bars; Table [2] and 



Barman et al.l (J2011al )) with atmospheric models (green, pink and blue curves) and predicted 
in-band fluxes for each photometry point (green, pink and blue horizontal lines, which span 
the filters' half- max profiles). The 850 K e quilibrium chemistry r nodel (green) is the best-fit 
"thick-cloud" model for HR 8799 b from iMadhusudhan et al.l (120 111 ). Two other models 
(pink and blue) suppress the CH4 mixing ratios and enhance the CO mixing ratios by 10 x 
and 100 X with respect to the equilibrium chemistry model. The H and K spectroscopy are 
well-fit by the non-equilibrium chemistry models. Our new 3.3yum photometry is best fit by 
the 100 X CO, 0.01 XCH4 model. However, this model predicts a fiux that is substantially 
higher than published measurements at 3.8/im (L'). None of the models are able to reproduce 
the relatively fiat SED from 3.3-3.8/im. 
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Fig. 7. — Same as Figure Elb ut with a mixed cloud mode l (93% 700 K A-type clouds and 7% 
1400 K AE-type clouds from iMadhusudhan et al.l (l201ll )). The fit adequately reproduces all 
photometry (except in the M-band filter, where additional CO absorption would rectify the 
discrepancy), and greatly fiattens the SED from 3.3-3.8/im compared to the non-equilibrium 
models shown in Figure |6l At H and K, the model generally reproduces the shape of HR 
8799 b's observed spectrum (i.e. negligible CH4 absorption), but is fiatter. 
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Fig. 8. — Same as Figure |6] (top) but for HR 8799 c and d. As was found for HR 8799 b, our 
non-equilbrium chemistry models are unable to fit the 3.3/xm-L' colors of HR 8799 c and d. 
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Fig. 9. — Same as Figure [7] (top) but for HR 8799 c and d. We purposely choose the same 
cloudy atmospheres to mix as were used for HR 8799 b. We are able to fit HR 8799 c and d 
with mixed cloud atmospheres, but using a higher mixing fraction of the 1400 K atmosphere 
than was found for HR 8799 b. The addition of non-equilibrium chemistry to these models 
would likely improve the fit. 
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Fig. 10. — Thick cloud atmosphere for HR 8799 e (plotted with the s ame symbols used in 
Figure s EMI). In this fit, we ignore the 3.3/im photometry as was done by lMadhusudhan et al. 
(|201lh when fitting HR 8799 b, c and d. 
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Fig. 11. — Same as Figures |6] (top) and |8] but for HR 8799 c and d. As was found for HR 
8799 b, c and d, our non-equilbrium chemistry models are unable to fit the 3.3/im-L' colors 
of HR 8799 e. 
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Fig. 12. — Same as Figures [7] (top) and [9] but for HR 8799 e. The mixing fraction of 
the two atmospheres is similar to what we found for HR 8799 c and d. The addition 
of non-equihbrium chemistry to these models would slightly improve the fit, although the 
photometric error bars for HR 8799 e are large enough that this is not necessary. 
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